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The problem of electrons energy spread in the linear back Compton scattering process has been 
considered in the paper. The adjoint kinetic equation for the electrons passing through a laser flash 
has been obtained and used to get the equations for the mean energy and the variance of the energy 
distribution. The equations for the distribution moments have been obtained and solved by the 
iteration method in the continuous slowing down approximation with approximate description of 
the energy loss fluctuation. It has been shown that the variance of the energy distribution as a 
function of the electron - photon collisions number in laser flash has a maximum for the electron 
beam with small incident energy spread. The beam energy spread damping in the laser-electron 
storage ring has been studied. The data of approximate analytical calculations are compared with 
the results of the Monte Carlo simulation. 
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I. INTRODUCTION 

The interaction of a high energy electron beam with 
an intense laser flash is now considered in such projects 
as compact x-ray sources, 7 — 7 colliders, diagnostics of 
low emittance electron beams, laser cooling and others 
0. The laser flash intensity in these projects may be so 
high that an electron can undergo several successive col- 
lisions passing through a photon bunch & B. It means 
the number of scattered photons is larger than one. Tak- 
ing into account the discrete character of the Compton 
back-scattering process, one can see an analogy between 
the electron passage through a laser flash and through 
a condensed matter. The stochastic theory of particles 
penetration through a matter can be found in the book 
Q . The model of multiple bremsstrahlung process which 
leads to emission of a few photons by a single electron was 
developed in Here we consider the multiple Comp- 
ton back-scattering process using the same approach. For 
simplicity we use the approximation of an uniform pho- 
ton concentration in the laser flash: 
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where A is the total energy of the laser flash; u>q is the 
photon energy, and II are radius and length of the 
photon bunch. In other words our approximation corre- 
sponds to the short laser pulse (7l <C zr, zr being the 
Rayleigh length). 

The kinetic equations have been obtained for monoen- 
ergetic electrons and for electrons with incident energy 
spread. These equations are transformed into the equa- 
tions for the distributions moments. The method for an 
approximate solution of the equations has been derived 
for the mean energy and for the variance of the energy 



distributions. It has been shown that the energy distri- 
bution variance as a function of the target thickness is a 
curve with a maximum if the incident energy spread is 
small. 

It is shown that the equation for calculation of the 
variance can be used to study the variation of the energy 
spread in the laser-electron storage ring (LESR) where an 
electron beam repeatedly encounters the laser pulses. It 
has been shown that the variation of the energy spread in 
this case depends on the incident value of the variance. 
The energy spread damping takes place if the incident 
energy spread is larger than some limit value. 

The results of analytical calculations are compared 
with the data of Monte Carlo simulation. 



II. TRANSPORT EQUATIONS 

Let us consider electrons with an incident energy Eq 
mc 2 travelling through a bunch of photons with an energy 
luq. The process of an electron transport in a light target 
is a random process where both the collisions number 
along the path I and the energy loss in each collision are 
random quantities. Therefore, the energy of electron e 
passing a laser flash is random quantity too. 

The angular deflection of an electron in the Compton 

back-scattering process is small (~ — -r) and may be ne- 

mc 2 

glected (the straight-ahead approximation). In this ap- 
proximation the probability density function P m (e|eo,Z) 
(the subscript "m" means "monoenergetic") describing 
the energy distribution of electrons travelling a path I 
in the photon bunch obeys the adjoint balance equation 
(the Kolmogorov-Chapman equation) j|, 



P m (e\e , 1) = (1 - sE(E ))P m (s\e , 1 - s) 
£(w;e ) 
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where £(£o) and £(w;£o) are the total and differential 
macroscopic cross-sections of the Compton scattering, s 
is a small part of I, and a; is the energy of scattered 
photon (0 < lo < u max ), 



= £o 



1 + x 



is the maximum value of the scattered photon energy, 

^ _ 4^q£q 

(toc 2 ) 2 ' 



(the subscript "c" means "convolution"). 

Eq. (J2J) can be transformed into the equation for the 
distribution moments 



£ rn( £ o,i)= / e k P m (e\e ,l)de. 



The equation for e^(e , |(| is 

d 

— e^(eo, I) + S(eo)£m(eo, 



S(e ) = 2n L a(£ ) ■ 
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ct and — — are the total and differential cross-sections, til 
dio 

is the concentration of laser photons in the bunch. 

Note that E(e ) is the mean number of collisions of an 
electron per unit path length and T,(io;e ) is the mean 
number of an electron collisions with the energy loss in 
unit interval about lo per unit path length. 

The first term in the right side of Eq. ([!]) corresponds 
to the electrons, which passes the path s without colli- 
sions, and 1 — sS(eo) is corresponding probability. These 
electrons have to lose the energy £0 — £ along the rest 
path I — s. 

The second term corresponds to the electrons, which 
undergo the first scattering passing the path s, and 
sE(£o) is the corresponding probability. The energy 
of the electron after the first scattering equals So — lo, 
where lo is a random energy of the scattered photon, and 
E(w; £o)/S(£o) is the probability density function of lo. 
In the limit s — > Eq. ([!]) gives the adjoint integro- 
diffcrential equation for the function _P m (£|£o, I) || n|: 

d 

— P m (s\e ,l) + E(£ )P m (s\e ,l) 



E(w;e o )J , m(e|eo-w,l)dw = (2) 



with boundary condition 

P m {e\e J)\i=o = S(e - £ ), 

5(e — £0) being the Dirac (5-function. 

If the electrons incident on the light target have an en- 
ergy distribution Po(s), the energy distribution of elec- 
trons behind the bunch is the convolution of Po(£) and 
-Pm(e|eo,0 : 



The boundary condition for the moments is £^(£o,Z) 

is 

e^(e 0l l)\i =0 = eg. 

The equation for the moments of the distribution 
P c {e\l) describing transport of electrons with an incident 
energy spread can be obtained from Eq. (||): 

00 00 

i£(0= / e k P c (e\l)de= J P (e')eZ(e' J)de' . (5) 



It is seen that they are expressed in terms of the moments 
for the monoenergetic electrons and the incident energy 
distribution. 

If the energy loss of an electron in one collision is small, 
the integro-differential equation for the moments (^) can 
be transformed by the Taylor expansion of integrands: 



e<f n (e Q -u},l) w £^(£0,0 - w^-£^(£ ,0 

This gives the partial differential equation: 

d — d — 

—4, (£ ,0 + 0{s o )—s^(e o ,l) 



1 d 2 

27(£o)^«4(eo,0 = 0. 



(6) 



The quantities /3(£o), and 7(£o) in @ are the moments 
of the macroscopic differential cross-section: 



(3(e ) = J loY>(lo; e )du , 




Pc(e|0 = [ Po(e')P m (s\e',l)de' 



(3) 



7(£ ) = J lo E(u);e )du>. 
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Two first terms in Eq. correspond to the continu- 
ous slowing down approximation, whereas the third term 
gives approximate description of the energy loss fluctua- 
tions. 

The Eq. (|) written for k = 1, 2 can be transformed 
into the equation for the variance 

CmOo, = £mOo, ~ £m(>0, 0- 

The equation has a form || 

d did 2 

—a m (e ,0 + P{e Q )—a m (eo,l) - ^M^o^feo, 



= 7(eo) (^^( e °'Q J -( ? ) 

The Eq. (|) for the energy distribution moments in the 
case of the electron beam with an incident energy spread 
can be simplified if the distribution P (e) is a symmetric 
function with maximum at a point Eq and small variance 



lu e 
where y = — and ro = - is the classical radius of 

€q TfLC 

electron |2j. 

In the energy region of our interest the invariant di- 
mensionless parameter x is small (x <C 1) and the inter- 
action coefficients fl(eo), and 7(20) are described by the 
approximate formulas 



7(eo) 
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-s x 



(11) 



(12) 



ttTq being the Thomson 



where So = ZnL&T, - 
cross-section. 

The quantities /3(sq) and 7(eo) are the mean energy 
loss and the mean squared energy loss of an electron per 
unit path length. 



(e - e ) 2 P (s)de. 



In this case the Taylor expansion of the function £ m (£, I) 
in d) makes it possible to express the moments of the 
distribution P c (e\l) in terms of the moments e^(£o,0 
corresponding to monoenergetic electrons: 



-.'■(/)= t(--<../> + y^4^o,0- (8) 



The equations (||) written for k — 1, 2 can be easily trans- 
formed into the equation for the variance 

a c 2 (o=i!(/)-^ 2 (0- 

The equation is 

+ y (^^( £o ,0 + 2(^^( £ o,0) 2 ) -(9) 

It describes the variance of the energy distribution for 
the electron beam with an incident variance uo in terms 
of stochastic characteristics for monoenergetic electrons. 



IV. SOLUTION OF EQUATIONS FOR 
MONOENERGETIC ELECTRONS 

The partial differential equations (||), (0) with the in- 
teraction coefficients ( |TT| ) , ( |l2| ) describe the mean energy 
and the variance of the energy distribution for inciden- 
tally monoenergetic electrons passing a path /. The equa- 
tions can be solved by the iteration method. In the first 
approximation, where the terms with the second deriva- 
tives are neglected, the solutions are 



i( £ 0,0 = 



£0 



1 + \nx 



Cm (£0,0 



7 e\x 2 n 
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(13) 



(14) 



where n = "EqI is the mean number of collisions in the 
light target H. 



The mean value (13) is a monotonically decreasing 
function of variable I. The Taylor expansions show that 
for small I (nx <C 1) 



III. COMPTON SCATTERING 
CROSS-SECTIONS AND INTERACTION 
COEFFICIENTS 

In this paper we restrict our consideration to the linear 
Compton scattering. The differential cross-section of this 
process for relativistic electron is 



da(y;x) 2?rr§ 



dy 

(i-y- 



4y 



1-y x(l-y) x 2 {l-y) 2 



,(10) 



£m(£o,0 ~ £o(l - -y) 



and for large I (nx » 1) 



e m {s ,l) w — 1- — 
nx \ nx 



Similarly, it can be shown that the variance (^) grows 
with I for thin targets: 



CT m( £ o,0 ~ —nelx 2 {\ - 2nx) 
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and decreases for large I: 

o , 28 £q . 8 

5 nrar ncc 

It follows from ( fTi]) that the variance af n (eo,l) has a 
maximum at the point, where 

2 1 
n = ■x- 

3 x 



or 



£m(£0) 



£0 



The energy spread damping for large I is due to the 
energy loss (3(e) per unit path length is proportional to 
e 2 , therefore, higher energy electrons in a bunch lose more 
energy than lower energy electrons (§, |9) . 

Using Eqs. (|J) and @ one can derive the formula 

o"m( £ o,0 = 77^— ( 1 - — ) , 
10 e V £ o/ 

which was earlier obtained in 0], 



V. SOLUTION OF EQUATIONS FOR 
ELECTRON BEAM WITH INCIDENT ENERGY 
SPREAD 

Substitution of Eqs. ©, @ in Eqs. (|) and §) 
makes it possible to calculate the mean energy and the 
variance of the energy distribution for the beam with the 
incident energy spread a$. 

These calculations show that the mean value eZ{l) is 
decreased monotonically with /: 

BS(0«*o f 1 ~?( 1 + 4) 
V z £ o 

for small I and 
for large Z. 

But the dependence of variance a 2 (I) on I is determined 
by the incident energy spread. The Taylor expansion of 
Eq. (ft) in powers of / gives 
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whereas for large / it is decreased with I: 



a 2 (e ,l) 
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It means that the variance (t 2 (0 has a maximum at some 



/ if a 2 < o 2 maxl where 



L 2 i 

40 X£ °l-§x' 



(15) 



For (To = c maa; this maximum is at the point 1 = 0. If 
do > (T ma2; the variance (Tc(0 i s a decreasing function of 



VI. MULTIPLE CROSSING OF LASER BUNCH 

It was pointed out in M, [l(| that effective radiative 
cooling of an electron beam can be realized in LESR 
where electrons repeatedly interact with an intense laser 
pulses. Eq. (||) can be used to study the variation of the 
energy spread for electrons multiply crossing the light 
target. If the lost energy is restored by an rf accelerating 
system after each turn, the variance of energy distribu- 
tion after n-th turn can be calculated using Eq. (^|) with 
<To replaced by the variance after (n-l)-th turn: 



where 



a 2 (n) 



1 / 8 2 



,(£0,0 +V v 2 (n- 1), 



77 = H^ <Tm(£o,0+2 W m(eo ' 0) ' 

It follows from (O) that for \r]\ < 1 



a 2 (oo) = a^(e , 0(1 + V + V 2 + V 3 + •••) 



(16) 



(17) 
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Let us point out that the multiple crossing of the laser 
bunch leads to the energy spread damping if the incident 
energy spread a 2 > a 2 (oo). In the opposite case (a 2 < 
a 2 (oo)) the energy spread is increased. 

Limit value of energy spread a 2 (oo) depends on the 
thickness of laser bunch especially for high energy elec- 
trons. For small I 



<r 2 (oo) « a 



2 

max 



, _ (15 + 282) 

1 — n x -, r 

(20- 21x) 



whereas for large I 



(T 2 (00) 
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It is easy to see that for small x 



<t(oo) 
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where Xl is the laser radiation wavelength and A e = — 

mc 

is the Compton wavelength of electron. This result coin- 
cides with that one given in p0| . But our approach allows 
to estimate the number of turns in a LESR to reach the 
limit value of the energy spread. 



VII. NUMERICAL RESULTS 

Numerical calculation of the mean energy and the vari- 
ance of the energy distribution were made for 100 MeV 
electrons interacting with 1.24 eV laser photons and for 
5 GeV electrons interacting with 2.48 eV photons. Eqs. 
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(||) and (0) were used to calculate the relative energy 

spread ^1,} ■ The results are given in Fig. |l| for 3 values 

of the incident energy spread. They agree with the data 
obtained using Eq. (15) of 11 



It is seen from the figure that the relative energy spread 
of the electron beam as a function of the collisions num- 
ber n is a decreasing function if the incident energy 
spread is large enough. But for the beams with small 
incident spread the "heating" of the beam takes place 
instead of " cooling" for small I . 

The dependence of the energy distribution variance on 
the number of turns in LESR is shown in Fig. ^ for two 
incident energy spreads. The calculations were made for 
LESR suggested for laser beam cooling in 10 . In this 



case wo = 1-24 eV, Eo — 100 MeV, the mean energy 
of scattered photon equals U — 2 £ ° x = ^ keV, the 
average ring radius equals 1 m, and the energy loss per 
turn Auj — 25 keV. The mean number of collisions in 
laser bunch n — 4^ = 0.27. It is seen from the figure that 
the asymptotic value of the energy spread is established 
in ~ 5 • 10 3 turns and it takes ~ 100 /isec. 

It is seen from Fig. || that the energy spread of elec- 
trons multiply crossing the light target is decreased or 
increased depending on its incident value. 

The limit value of the relative energy spread (j(oo)/eq 
as a function of the light target thickness is given in Fig. 
||. It is seen from the figure that the variance of the 
energy distribution is decreased with the target thickness 
increasing. 



VIII. MONTE CARLO SIMULATION 

It is supposed in our Monte Carlo simulation of the 
back Compton scattering that the incident energy distri- 



bution of electrons is Gaussian with given variance Cq. 
The number of electron collisions with laser photons is 
supposed to be random. It is selected from the Poisson 
distribution with fixed n. The simulation of individual 
collisions is carried out in the electron rest frame using 
the Klein and Nishina formula with the Lorentz transfor- 
mation to the lab system. In the same way as in analyti- 
cal calculations above we neglected the angular deflection 
of electrons but accounted for the energy decreasing after 
each collision. 

Fig. |^ shows the energy distributions of 5 GeV elec- 
trons for several values of the light target thickness. It 
should be pointed out the discontinuous of the spectra for 
small n due to the single scattered electrons. It is seen 
from the figure that the energy spread increases with the 
target thickness increasing for small n and decreases for 
such n, which are greater than determined by Eq. 

©• 

Fig. |^ and || shows the electron energy distributions in 
LESR for several number of turns. 

All results are obtained with statistics more than 10 6 
trajectories. 

It should be pointed out that the analytical formulas 
obtained in the continuous slowing down approximation 
with approximate consideration of the energy loss fluc- 
tuations are inapplicable for highly relativistic electrons 
if the mean energy loss in one collision is comparable to 
Eq . The Monte Carlo method has to be used in this case. 

Fig. shows the relative energy spread as a function of 
the light target thickness calculated in continuous slowing 
down approximation and by the Monte Carlo method for 
250 GeV electrons. Fig. || shows the energy spectra of 
these electrons. Notice that the Monte Carlo technique 
allows to get the results with taking into account the 
energy dependence of the interaction cross-sections and 
the lateral distribution of photons in the laser flash. 
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FIG. 1: The relative energy spread a c (l)/e c {l) as a function of the target thikness for various initial energy spread. 
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FIG. 2: The relative energy spread as a function of the number of turns in LESR, £o = 100 MeV, uo = 1.24 eV, and n = 0.27 
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FIG. 5: Evolution of the energy distributions P(e) for electrons with ao < <r(oo) multiply crossing the light target. Solid line 
is the Gaussian distribution with the variance cr(oo). Number of turns is shown near the curves. 




FIG. 6: Evolution of the energy distributions P(e) for electrons with <to > cr(oo) multiply crossing the light target. Solid line 
is the Gaussian distribution with the variance <r(oo). Number of turns is shown near the curves. 
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FIG. 7: The relative energy spread c c (/)/ej(Z) as a function of the target thickness. Solid line - Eqs. ([]), (|^), points - Monte 
Carlo simulation. 




FIG. 8: The energy distributions of initially monoenergetic electrons for n = 



0.1, 0.5, 1 and 2. 



